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The van der Wahs interaction between a carbon nanotubes sheet (CNTS) and a rare gas atom, is 
studied using both atomistic and continuum approaches. We present analytical expressions for the 
van der Waals energy of continuous nanotubes interacting with a rare gas atom. It is found that 
the continuum approach dose not properly treat the effect of atomistic configurations on the energy 
surfaces. The energy barriers are small as compared to the thermal energy, which implies the free 
motion above the CNTS in heights about one nanometer. In contrast to the energy surface of a 
graphene sheet, the honeycomb lattice structure in the energy surface of a CNTS is imperceivable. 
Defects alter the energy surface which therefore influence the gas absorbtion mechanism. 



I. INTRODUCTION 

In the past two decades, various properties of carbon 
nanotubes such as their unique structure and mechan- 
ical, chemical and electronic properties, have been at- 
tracted much attention [1]. Zhang et al used a solid 
state sheet fabrication technique and could synthesize 
highly oriented free standing carbon nanotubes sheets 
(CNTS) [2]. The produced CNTS combines high trans- 
parency with high electronic conductivity and are highly 
flexible and provide giant gravimetric surface areas. The 
long, high perfection nanotubes are needed to maximize 
the electrical and thermal conductivities and mechanical 
properties [2] . Recently Li et al [3] have fabricated highly 
engineered ultra-thin anisotropic single wall carbon nan- 
otube films formed on three dimensional micro-patterned 
substrates. Therefore CNTS can be considered as a real 
multifunction sheet for the future technologies. 

The two-dimensional structure of CNTS makes it pos- 
sible to examine the possibility of the two-dimensional 
motion of an external particle over them [4, 5]. The ab- 
sorbtion of rare gas on the CNTS can change the lat- 
tice parameter of nanotubes. In this regard Bienfait et al 
found experimentally the dilation of the apparent carbon 
nanotube bundle in the presence of argon gas atoms 

The potential energy surface (ES) of CNTS explains 
different motion-related phenomena at nanoscale as well 
as the various directed motions on the carbon nanotube- 
based motors jgl, Since rare gas atoms are experi- 
mentally used to control temperature, pressure and other 
thermodynamical properties of nanoscale materials [l, [9| , 
it is important to investigate the van der Waals (vdW) 
potential energy surface around a CNTS. Recently, we 
found a diffusive motion for the motion of Ceo molecule 
over graphene sheet and depicted that the energy barri- 
ers are small as compared to the thermal energy [6[. By 
employing a different approach, i.e. Langevin dynam- 
ics, we showed that a rare gas over graphene sheet has a 
diffusive motion [foj . 

There are two common approaches for studying the 
various properties of carbon nanotubes based nanostruc- 
tures: 1) atomistic and 2) continuum. Yakobson et al 
[ill have shown that the mechanics of carbon nanotubes 



is interpretable using the continuum elasticity of shells. 
The pioneering work of Yakobson and co-workers made 
common practice of using of continuum models for sim- 
ulations of mechanical properties of nanotubes jl2|-[l3- 
In recent studies, it has been shown that atomistic and 
continuum approaches can be combined to geth er to study 



mechanical behavior of different materials [12|. However, 
in the present study, we will show that employing of the 
continuum approach for energy based properties (but not 
elastic energy) is not appropriate. 

In this study the ES of a CNTS due to the vdW in- 
teraction with a rare gas atom is investigated. We show 
that the ESs are periodic in the planes which are paral- 
lel to the CNTS layers. We find that there are signifi- 
cant differences between ES of the atomistic model and 
those predicted from the continuum model. The binding 
energies and the physical-bond distances obtained from 
the atomistic model differ from the continuum model. 
We conclude that the continuum model is inappropri- 
ate approach for energetic consideration and yields inac- 
curate results for the chemical absorbtion of atoms into 
the CNTS. However mechanics of atomistic systems (car- 
bon nanotubes) greatly benefits from continuum mechan- 
ics [nl, [l2|. In addition, we have investigated the ES of 
graphene and compared it to the one obtained for CNTS. 
It is found that the honeycomb lattice feature, which is 
observed in the ES of graphene, does not appear in the 
ESs of CNTS. This is attributed to the large curvature 
of carbon nanotubes in the CNTS sheet 

Results for the CNTS comprises the zigzag nanotubes 
differ from the CNTS made of other types of nanotube 
chiralities. From our barrier energy calculations, it is 
found that the barriers for the motion of the rare gas 
atom on a CNTS become shallow at positions away from 
the CNTS. Moreover, the components of the force are 
calculated and are compared to those obtained from the 
continuum model. We show that the simplest kind of 
defect (i.e. vacancy), destroys the periodicity of the ES 
and reduces the energy barriers. 

This paper is organized as follows. In Sec. II and HI 
the details of atomistic and continuum models are pre- 
sented Sec. IV contains main results including those for 
the energy surfaces, force calculation, numerical solution 
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for the equation of motion and the results for ES of the 
defected CNTS. In Sec. V the results are summarized. 



II. THE ATOMISTIC MODEL 

A CNTS comprises many layers which are stacked on 
top of each other and are along z-axis. Many aligned 
carbon nanotubes which are parallel to the y-ax.is con- 
struct each layer. The carbon-carbon bond length, ao, 
IS 1.42 A. In our model, the origin of the xyz Cartesian 
coordinate system is set to be at the center of the top 
layer. Carbon nanotubes form different types, which can 
be described by the chiral vector with two indexes (m, n), 
where m and n are integers of the chiral vector equation 
^ = m~^i + n~^2- Here and are two basic vec- 
tors of graphene lattice and are given by "^i = \^aoi 
and = A/3/2ao^ + 3/2aoj, respectively. If m = 0, 
the nanotubes are called zigzag (ZZ). If n = m, the 
nanotubes are called armchair (AC). Otherwise, they are 
called chiral (CR). The radius of a carbon nanotube with 

(m, n) indexes is given by Rm,n = ^^^^-y/rri^ + + mn. 

Figure [T]^a) depicts eight aligned (5,5) carbon nan- 
otubes from the top layer. Graphene is a planar struc- 
ture made of carbon with one-atom thickness, where 
each carbon atoms have sp^ hybridization with their 
neighboring atoms and are densely packed in a honey- 
comb crystal lattice. Carbon nanotubes can be formed 
by rolling up graphene sheets. In order to study the 
van der Waals interaction between an argon atom and 
a CNTS, we have employed the Lennard- Jones (LJ) po- 
tential. The LJ potential provides both the repulsive 
and attractive nature of the interaction between an un- 
charged noble gas and each atom of the CNTS. The LJ 
potential is the widely used potential in various simu- 
lations [15j for two interacting uncharged particles, i.e. 
u(t) = ^eiialvY^ — {a/r)^)^ where r is the distance 
between two atoms, e is the depth of the potential well, 
and <j is the distance at which the potential becomes zero. 
To model the interaction between two different types of 
atoms such as carbon and argon, we adjust LJ param- 
eters using the equations e = and <j = ^^^^^ . 
We use the parameters of carbon and argon atoms as 
(Ji = 3.369 A, Ex = 2.63 meV, G2 = 3.405 A and 
82 = 10.23 meV, respectively. 

Total potential energy stored between an argon atom 
and a CNTS is written as the sum over all contributions 
of the carbon nanotubes in each layer. 



N/2 



K M 



i = -Ar/2t=0j=l^=l 



12/1 



(1) 



where 8r = 1 1^ — (1^* ^ —tD k)\ and 1^ = x y j -\- zk 



refer to the position of the argon atom and 1^* ^ is the 



position vector of j^^ carbon atom in the i^^ carbon nan- 
otube (of the t^^ layer). Each nanotube has M atoms 
and the total number of carbon nanotubes in each layer 



is + 1. The sum over / also, is written for varying 
of the powers between 12 or 6. The sum over t counts 
the number of layers {K is the total number of layers) in 
different height along z-axis, such that each layer is sep- 
arated by D from its neighboring layer. Note that each 
layer is shifted along the x-axis with respect to its top 



and bottom neighbors, i.e. r 



± |z where 



d is the lateral distance between tubes in each layer, see 
Fig. [H In the CNTS having many layers, we show that 
the main contribution of the energy comes from the top 
layer, thus setting K = 1 would be a good approxima- 
tion. We will prove this effect in Sec. IV-A (henceforth 
we set i^T = 1). Note that the potential energy in Eq. ([T]) 
is periodic along x-axis at fix z independent of the chi- 
rality and z, but the periodicity along y-ax.is depends on 
the chirality, i.e. Erix^y^z) = Et{x /j^d^y -\- u Sm,n^ z) 
where fi and u are two integer numbers and Sm,n is the 
periodicity number along ?/-axis, e.g. ^5^5 = ^/Sao and 
(^9 = 3ao. In fact, the periodicity of the ES is a conse- 
quence of the crystalline structure and parallel alignment 
of the carbon nanotubes. Moreover, the gas absorbtion 
into the CNTS depends on the ES periodic structure. 

The gradient of the potential is the force experienced 
by the argon atom 
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Newton's second law should be solved numerically for 
determining the path of the motion. 



III. THE CONTINUUM MODEL 

Considering a hollow and continuous cylinder as a 
simple model for the carbon nanotube is the main 
assumption for the continuum approach. Again, A^+1 
aligned tubes along ?/-axis are separated by the distance 
d in each layer of CNTS. The obvious difference between 
the atomistic model and the continuum model is the 
absence of the chirality effect in the continuum model. 
The other important difference is the periodicity of ES 
along ?/-axis in the atomistic model (not in the continuum 
model) which will be discussed further in the next sec- 
tion. The LJ potential energy for A^+1 cylinders which 
are interacting with a rare gas atom and are parallel and 
aligned along the ?/-axis (when K = 1) is 
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Ec{x,y,z) = Ei{x,y,z), 

i=-N/2 
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where Ei{x,y^z) is the contribution of the i^^ nanotube 
and is given by 



2 /.27r rL 



Ei{x,y,z) = E / / 
i^i Jo Jo 



2nii-lY+' R„,md(/>dyi 
-2S'm(<^i) + (2/-2/i)2]6/i' 



(4) 



3 



m ik M ^ ilk M ^ 




(x,z) 




(a) 



(b) 



FIG. 1: (Color online) (a) The atomistic model of a CNTS. The top layer of an AC-CNTS comprises eight nanotubes with 
M = 400 and (i=10A. (b) The continuum model of a CNTS. The front view of three continuous tubes with the relevant 
variables for Eq. (|4j). The red sphere is an argon atom located at {x,y, z) position. 



with ni = 4£:cr^^/^I] where S is the mean surface density 
of carbon atoms for a graphene lattice which is S = 
. In Eq. (jl]) {xi^yi^Zi) is the coordinates of the 

chosen element from the i^^ nanotube, L is an arbitrary 



length, i?^ 



xj + is the radius of each nanotube, 



tan-i(^), 



{x -\- id) -\- z . Note that "z" is a positive integer 
number where Xi < and vise versa. In Fig. [TJb) we 
illustrate the front view {in x — z plane) of only three 



nanotubes from the top layer. The relevant variables 
for the model are exhibited in Fig. [IJb). The following 
integrals are taken over the middle nanotube referring to 
z = l in Eq. (g]). 

The integrals over the longitudinal viable {yi for the i^^ 
cylinder) are done analytically and the integrals over the 
azimuthal angle {(j)) are obtained numerically by employ- 
ing Mathematica software. The result of the integration 
over yi for repulsive term is 
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and for attractive term is 
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where A = — 2 Sin{(pi) and C = V ~ Vi- Combining two above integrals, the continuum model energy is written as 
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Differentiating the above equation gives the force ex- perienced by the argon atom due the interaction with i 



th 
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tube: 

where Ar = {x -\- id — Xi)i -\- {y — yi)j -\- {z — Zi)k. 

We will compare the results of the continuum model 
with those of the atomistic model in the next section. 



IV. RESULTS AND DISCUSSIONS 
A. Energy surface 

Figure [2] shows two energy surfaces (ES) in the x — y 
plane above the CNTS at z = 6.5 A which are related 
to the atomistic (Eq. ([1])) and the continuum (Eq. (|3|)) 
models, i.e. Erix^y^Q.b A) and Ec{x^y^6.5 A) respec- 
tively. The ESs were calculated with K=l^ N = 7 and 
d=10 A (these are fixed in this study). For the atomistic 
model (Fig.[2](a)) we set M = 400 for AC nanotubes with 
(5, 5) indexes. The colored tubes indicate the location of 
the top layer of the CNTS (includes eight aligned carbon 
nanotubes). We found that the chosen values for M and 

are sufficiently large to eliminate boundaries when ES 
is considered in the range of \x^y\ < 20 A. In contrast 
to the ES obtained for the graphene [6], from the ES 
in Fig. [2fa) one can not identifies the honeycomb lattice 
structure of the nanotubes. This is because of the curva- 
ture of the nanotubes. Moreover, the ES approaches to 
a fiat surface far away from the sheet. 

As can be seen from Fig. [21 the main difference between 
two models is the appearance of the periodic peaks along 
7/-axis in the atomistic model. In spite of the contin- 
uum energy surface, i.e. Ec^ Et has many peaks/ valleys 
along the 7/-axis (above (id^y) points for each integer 
'i'). These valleys make the connection between two ad- 
jacent channels easier and will also enhance the proba- 
bility of the lateral motion (along x-axis) of the particle 
through the mentioned valleys. Furthermore the appear- 
ance of the peaks along the x = id/2 line can be un- 
derstood regarding to the closed energy curves around a 
single nanotube (Fig. [H^a)) in the x — z plane. Assume 
two such nanotubes get laterally closer from infinite dis- 
tances. At certain distances due to the overlap of their 
energy curves, the mentioned peaks appear. 

We depict ES of ZZ-CNTS with (9,0) indexes and CR- 
CNTS with (6,4) indexes in Fig. [3] (we used shorthand 
letters ZZ, AC and CR for ZZ-CNTS, AC-CNTS and CR- 
CNTS, respectively in the figures legends). Note that 
Rg^o ^ i?6,4 — ^5,5 = 3.39 A and that is the reason for 
choosing (9,0) and (6,4) indexes. It is clear from Fig. [3] 
that the periodicity of Et along y-ax.is depends on the 
(m^n) indexes. 

Increasing the number of layers, i.e. varies ES 
slightly implying that the major part of the ES is due 
to the interaction with the top layer. In order to prove 
this fact, we show the variation of Et for AC-CNTS 



with K= 3 and K= 1 versus z in Fig. |4j The ratio 
Et\k = 3/ Et\k = 1 varies around one, thus setting K = l 
is a reasonable assumption. The same behaviors were 
obtained for the variations of the mentioned ratio versus 
X and y. 

Figure [5] shows the variation of ES for an AC-CNTS 
along z-axis at different (x, y) points, namely £^t(0, 0, z), 
ET{d/2,0,z) and Ec{0,0,z), Ec{d/2,0,z). Although 
these curves are similar to u{r) function, but indeed 
they are different. The curves for ZZ-CNTS have dif- 
ferent minima as compared to those for AC-CNTS and 
CR-CNTS. The minimum of energy curve for ZZ-CNTS 
(above {x = O^y = 0) point) is around 6.99 A — Rq^q 
and is about 6.75 A-i?5,5 for AC-CNTS and CR-CNTS. 
The lowest energy occurs for Ec- In general, the en- 
ergy curves above {id/2^y) points have a shift with re- 
spect to those above {id^y)^ see Fig. [5l The appearance 
of the minima above {id/2^y) points are apparently in 
contrast to the appearance of the ES's peaks above the 
same points in Fig. [2] and Fig. [3l In fact, this is an 
unpredictable effect that when z decreases, the energy 
increases above {i y) and decreases above the points 
between the tubes (i.e. {id/2^y)) independent of the y 
values. Above {id/2^y) points the minima are due to 
the free space between adjacent tubes, i.e. the grooves 
between two adjacent nanotubes. Therefore the particles 
will prefer to aggregate between the tubes than above the 
tubes. Notice that above {i d/2^ 0) points the binding en- 
ergy (minima of the curves) is almost twice smaller than 
the one above (i(i, 0). The minima of the curves above 
(0,0) are around 70meV and above {id/2^0) points are 
about 128 meV which are 2.7 and 4.9 times larger than 
the thermal energy at room temperature, respectively. 
The binding energy of an argon atom above a graphene 
sheet is around ^ 85 meV [10| , which is smaller than the 
binding energy of an argon atom located at {id/2^0) 
above a CNTS. 

To investigate the periodicity of energy along x and y- 
axis the variations of ES are shown in FigEl Notice that 
the variations of energy versus x are almost the same for 
AC, CR and Ec while they are different from ZZ. In all 
cases the minima of the curves are equal, independent of 
the chirality an z values. Figure [6l^b) shows the variations 
of energy versus at z = 6.5 A. In spite of the curves in 
Fig.[6fa), the periodicity along y depends on the chirality. 
The Ec is a constant along y-axis. 

The maximum barriers are found for ZZ-CNTS and 
the minimum are for AC-CNTS. The barriers in the CR- 
CNTS vary with y resulting the possibility of the absorb- 
tion of external atoms in different energy levels. However, 
the barriers in the AC-CNTS are constant for fixed 
which are consequence of the sine shape of the curves. 
Figure [6fc) shows the variation of the energy barriers 
{5 = Emax — Emin) vcrsus z in log- log scales. The bar- 
riers decrease non-linearly with z which implies that far 
from the sheet, the particle has free-like motion. 
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FIG. 2: (Color online) (a) The vdW energy surface of a CNTS comprises of aligned armchair nanotubes interacting with an 
argon atom (located at 2; = 6.5 A), (b) The vdW energy surface of continuous tubes interacting with an argon atom (located 
at z = 6.5 A). 




(a) 



(b) 



FIG. 3: (Color online) The vdW energy surface of chiral nanotubes with (6,4) indexes in (a) and zigzag nanotubes with (9,0) 
indexes in (b), which are interacting with an argon atom (located at Z — 6.5A). 



B. Force calculations 



Using Eq. (|2f and Eq. ^ we calculated the various 
components of the force experienced by an atom above 
CNTSs with different chiralities. Here we report the 
variation of Fz versus Fx versus x and Fy versus y. 
Figure [Tl^a) shows the variation of the Fz versus z at 
two points (0,0) and ((i/2,0). Clearly the ZZ-CNTS at- 
tracts particle differently. The forces on the argon atom 
at (i d/2ft) point is more attractive than the one at {i dfi) 
implying tendency of the particles to aggregate between 
the tubes. 



Figure El^b) shows the variation of the F^ versus x for 
various 7/ at z = 6. 5 A. As we see the curves for the 
AC-CNTS and the CR-CNTS are similar. Figure [T^c) 
shows the variation of the Fy for an AC-CNTS versus y 
for different x and z = 6. 5 A. The absolute values of the 
forces decrease with increasing x. 



C. Equations of motion 

In order to find the path of the motion of an ar- 
gon atom above the CNTS, we solved equations of mo- 
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FIG. 4: (Color online) Comparison of the energy curves (as 
a function of x) produced by three layers of AC-CNTS (solid 
curve) and by one layer of AC-CNTS (dashed curve) for ^ = 
and z = 6.5 A. 




z(A) 

FIG. 5: (Color online) Variation of the energy versus z for 
different type of carbon nanotube sheets above two points 
(0, 0) and {i d/2, 0) in the x — y plane. 



tion numerically. The aim is to find the influence of 
the ES on the path of motion of an external atom. 
The study of an argon atom close to a single nan- 
otube surface provides further insights into the prob- 
lem. Therefore it would be helpful to show how ES en- 
ergy curves varies around an individual carbon nanotube. 
Figure [Hl^a) shows the energy curves for three z values, 
I.e. z= 6.25 A, 6.5 A, 6.75 A, 7 A. These cycloidal closed 
curves have been obtained by transforming the Carte- 
sian coordinates (in Eq. ([T])) to the polar coordinates. It 
is observed that with increasing the curves close to 



the circular shape which results in the flattening of ES 
at distances far away from the sheet 

The numerical solutions is presented for both the mo- 
tion of an argon atom around an individual AC carbon 
nanotube and an AC-CNT. The velocity Verlet algorithm 
is used for the numerical integration and the time is step 
set to be 1 fs. To solve the Newoton's equation one needs 
two initial conditions, namely the initial position of the 
argon atom and its initial velocity. We set the initial 
position to (0,0,6.5 A) and the velocity vector to an ar- 
bitrary vector ^ = (vx^Vy^Vz). We found that the larger 
initial velocities (corresponds to the higher kinetic en- 
ergy of the atom) give the larger probability for passing 
the barriers. Figure [8l^b) shows the front view of the two 
typical paths with ^ = (3,3,3) and ~^ = (5,5,5). It is 
also found that argon atoms with larger initial velocity, 
have larger fluctuation amplitude around the tubes. 

In addition, we give the results for the path of the mo- 
tion of an argon atom above an AC-CNTS. Figure [Sfc) 
shows the front view of the resultant. Notice that the 
particle prefer to aggregate within the grooves. Increas- 
ing the initial velocity of the particle enables the particle 
to overcome the barriers and visit more points at higher 
heights. In Fig. [8fb) and Figure [Sfc) the dashed inner 
circles indicate the carbon nanotubes locations. 



D. Defected CNTS 

Removing an atom from a perfect carbon nanotube 
makes CNTS defective. The ES varies in the presence of 
the defect. Defects changed the absorbtion mechanism 
and the paths of the motion of external atoms. Figure [9] 
shows the ES of an AC-CNTS (with (5,5) indexes) where 
vacant site was put in the middle of the central nanotube. 
Notice that the periodicity of the ES is altered around the 
vacancy. Vertical yellow arrow refers to the vacant site. 
We show the variation of Et{x^ z) versus both x and y 
in Fig.[9l^b). The solid curve corresponds to the variation 
versus x and the dashed curve is related to the variation 
versus y. Comparing these curves with those shown in 
Fig.[6fb,c) for perfect CNTS implies that the energy bar- 
riers are reduced due to the presence of the defect (gives 
more free motion along the x-direction). We have per- 
formed the same calculations for the CNTS with different 
chirality and found similar behaviors. Moreover, solving 
equations of motion in this case showed that the regu- 
lar aggregation between and above the tubes is missed. 
Increasing the number of vacant sites causes aperiodic 
energy surfaces. 



V. CONCLUSIONS 

Both the atomistic and the continuum models were em- 
ployed to study the van der Waals energy surface around 
a carbon nanotube sheet. The continuum model can not 
reveal the atomistic features of the problem. The car- 
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FIG. 6: (Color online) The periodic vdW potential energy versus x for (x, 0,6. 5 A) in (a) and versus y for (0,2/,6.5A) in (b). (c) 
Variation of the energy barriers (6 = Emax{x,y, z) — Emin{x,y, z)) along z. 




FIG. 7: (Color online) (a) Variation of the 2;-component of the force versus z which is acting on an argon atom located above 
(0,0) point, (b) Variation of the x-component of the force versus x which is acting on an argon atom located at (x, 0, 6.5 A), (c) 
Variation of the ^/-component of the force versus y which is acting on an argon atom located at (0,?/,6.5A), (1, 2/, 6.5 A) and 
(3,^,6.5 A). 




FIG. 8: (Color online) (a) Front view of the energy curves at various heights (b) and the path of the motion of an argon atom 
(with two different initial velocities) around an armchair nanotube. (c) The path of the motion of an argon atom with initial 
velocity (2,2,2) in semispace above a CNTS. 
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FIG. 9: (Color online) (a) Energy surface for an AC CNTS with a vacancy in the center of the middle tube. The yellow vertical 
arrow indicates the vacant site, (b) The variation of the energy versus x and respectively. For the dashed curve x = and 
for the solid curve y = for z = 6.5 A. 



bon nanotube sheet comprises zigzag carbon nanotubes, 
shows different periodicity in the energy surface as com- 
pared to those obtained from armchair and chiral car- 
bon nanotube sheets. The equation of motion was solved 
numerically and the path of the motion of a rare gas 
atom moving within the energy surface was determined. 



The continuum model predicts the constant energy along 
tubes axes which is in contrast to the results obtained 
from the atomistic model. The energy barriers decrease 
rapidly with height. Defects disturb the periodicity of 
the energy surface, reduce the barriers and change the 
gas absorbtion mechanism. 
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